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The  decay  and  annihilation  cross  sections  of  dark  matter  particles  may  depend  on  the  value  of 
a  chameleonic  scalar  field  that  both  evolves  cosmologically  and  takes  different  values  depending  on 
the  local  matter  density.  This  possibility  introduces  a  separation  between  the  physics  relevant  for 
freeze-out  and  that  responsible  for  dynamics  and  detection  in  the  late  universe.  We  investigate  how 
such  dark  sector  interactions  might  be  implemented  in  a  particle  physics  Lagrangian  and  consider 
how  current  and  upcoming  observations  and  experiments  bound  such  dark  matter  candidates.  A 
specific  simple  model  allows  for  an  increase  in  the  annihilation  cross  section  by  a  factor  of  106 
between  freeze-out  and  today,  while  more  complicated  models  should  also  allow  for  scattering  cross 
sections  near  the  astrophysical  bounds. 


I.  INTRODUCTION 


The  particle  physics  properties  of  dark  matter  are  important  for  three  distinct  aspects  of  its  behavior:  they  determine 
how  the  initial  abundance  of  dark  matter  arose,  they  govern  how  the  dark  matter  distribution  evolves  and  influences 
structure  formation,  and  they  delineate  the  possible  ways  in  which  dark  matter  may  be  detected.  Of  course,  these 
three  roles  are  not  typically  independent,  since  they  all  depend  on  the  prescribed  interactions  between  the  dark 
matter  particles  themselves  and  also  between  dark  matter  and  the  Standard  Model.  These  connections  often  provide 
a  powerful  motivation  for  particular  dark  matter  candidates  -  for  example,  the  freeze-out  abundance  of  weakly 
interacting  massive  particles  points  to  new  physics  at  the  weak  scale,  which  in  turn  leads  to  an  attractive  connection 
between  dark  matter  and  proposed  solutions  to  the  hierarchy  problem,  such  as  weak-scale  supersymmetry. 

The  idea  that  dark  matter  could  have  interactions  of  astrophysically  interesting  magnitude  has  received  a  good  amount 
of  attention  m,  motivated  in  part  by  purported  discrepancies  between  the  standard  ACDM  model  and  observations 
of  structure  on  small  scales  (as  described  in  [S] ,  for  example) .  While  most  approaches  of  this  form  concentrate  on  giving 
an  appreciable  scattering  cross-section  to  the  dark  matter,  it  is  also  interesting  to  consider  enhanced  annihilation  cross 
sections  [SJ. 

One  obstacle  to  simple  implementations  of  this  idea  is  that  the  required  cross  section  for  a  thermal  relic  to  obtain  the 
right  relic  abundance  is  close  to  the  weak  scale,  far  too  small  to  be  relevant  to  dynamics  in  the  late  universe.  In  this 
paper  we  explore  the  idea  that  the  dark  matter  cross  section  might  be  much  larger  now  than  it  was  at  freeze-out,  due 
to  the  evolution  of  a  background  field. 

In  a  cosmological  context,  the  evolution  of  background  fields  can  assert  a  significant  influence  on  the  properties  of 
dark  matter  as  a  function  of  spatial  location  or  cosmic  epoch  [TUHTO],  A  straightforward  way  to  achieve  such  effects 
is  to  invoke  a  light  scalar  field  that  interacts  with  dark  matter  and/or  ordinary  matter  as  well  as  through  its  own 
potential,  and  whose  expectation  value  feeds  into  the  dark-matter  properties.  A  popular  scenario  along  these  lines  is 
the  “chameleon  mechanism,”  which  acts  to  screen  light,  cosmologically  relevant  degrees  of  freedom  to  protect  them 
from  precision  local  tests  of  gravity  p2UH24| . 

In  this  paper  we  investigate  dark  matter  that  interacts  through  a  gauge  symmetry  with  a  coupling  constant  that 
depends  on  a  chameleonlike  scalar  field.  (The  effects  of  chameleon  vector  bosons  on  laboratory  experiments  were 
considered  in  [25].)  Just  as  the  properties  of  a  cosmologically  relevant  scalar  can  be  drastically  modified  in  the 
presence  of  local  density  inhomogeneities  or  after  evolving  over  cosmic  time,  so  the  interactions  of  dark  matter  may 
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be  modified.  We  are  able  to  find  a  model  in  which  the  late-time  interaction  strength  is  considerably  higher  than  that 
at  freeze-out  -  although  admittedly,  this  behavior  does  not  seem  generic. 

We  begin  by  reexamining  the  conventional  story  of  dark  matter  freeze-out  according  to  the  Boltzmann  equation, 
but  with  the  additional  ingredient  that  the  dark  matter  properties  are  evolving  with  time.  We  then  look  at  specific 
models  featuring  a  Dirac  dark  matter  particle  and  a  U(l)  gauge  symmetry  that  is  spontaneously  broken,  along  with 
a  chameleon  scalar  field.  We  study  the  cosmological  evolution  of  this  coupled  system  and  calculate  the  dark  matter 
properties,  including  annihilation  and  scattering  cross  sections.  Finally  we  exhibit  numerical  solutions  to  a  specific 
model,  showing  that  the  annihilation  cross  section  can  increase  substantially  during  cosmic  evolution. 


II.  THE  GENERAL  PICTURE:  EVOLVING  DARK  MATTER  IN  THE  EARLY  UNIVERSE 


Before  discussing  specific  models,  let  us  first  consider  how  the  usual  story  of  dark  matter  freeze-out  might  be  modified  if 
the  annihilation  cross  section  depends  on  the  dynamics  of  another  field.  In  the  next  section,  we  will  explore  Lagrangians 
that  couple  the  dark  matter  to  a  scalar  field  that  affects  its  interaction  cross  sections.  For  simplicity  we  work  in  a  flat 
Friedmann- Robertson- Walker  (FRW)  universe,  described  by  the  metric  ds2  =  — dt 2  +a2(t)  x  ( dx 2  +  dy 2  +  dz 2),  with 
scale  factor  a(t). 


The  decoupling  of  dark  matter  takes  place  in  the  early  universe  in  the  radiation-dominated  regime,  in  which  particles 
with  masses  m<T  are  the  dominant  component  of  the  cosmic  energy  budget.  To  a  good  approximation,  we  may 
therefore  ignore  contributions  from  nonrelativistic  species  in  thermal  equilibrium  with  the  radiation  and  approximate 
the  energy  density  as 


and  the  entropy  density  as 


where,  as  usual, 
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and  gi  is  the  number  of  internal  degrees  of  freedom  for  particle  species  i. 


(1) 

(2) 

(3) 

(4) 


For  T  >  300  GeV,  g*s  =  3*  =  106.75,  which  includes  all  particles  in  the  Standard  Model.  When  100  MeV  >  T  > 
1  MeV,  the  electron  and  positron  are  relativistic  and  so  g*s  =  g*  =  10.75.  At  the  temperature  of  the  CMB  today, 
T0  =  2.725  K,  g*s,o  =  3.91,  and  y*;o  =  3.36. 


Consider  a  dark  sector  that  was  in  thermal  equilibrium  with  the  visible  sector  at  some  very  high  temperature  scale, 
below  which  they  decouple  effectively  enough  to  consider  each  sector  separately  to  be  in  equilibrium.  The  visible 
sector  is  at  temperature  T  with  entropy  density  s(T),  while  the  dark  sector  is  at  temperature  Td  with  entropy  density 
sd(Td).  The  expansion  of  the  universe  is  governed  by  both  sectors  with 

gi°\T)=g*{T)  +  gi{Td)(^jA  ,  (5) 

but  quantities  in  the  dark  sector  (for  example,  the  dark  matter  annihilation  cross  section  and  number  density)  are 
determined  by  Td  [2fi]. 


Since  the  entropy  in  each  sector  is  conserved  independently,  the  assumption  that  the  two  sectors  were  in  equilibrium 
at  some  unification  scale  at  time  tu  allows  us  to  express  the  dark  bath  temperature  in  terms  of  the  visible  bath 
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temperature  at  some  later  time  t  via 


Td(t)  =  ats^n) 
g*s(t)T3(t)  g*s(tu) 


(6) 


All  Standard  Model  particles  contribute  at  tu  to  give  g*s(tu)  =  106.75,  and  all  dark  particles  contribute  to  <7*s(tu). 
In  what  follows,  we  will  use  the  temperature  of  the  visible  sector  and  convert  Td  to  T  as  needed.  For  convenience  we 
write 


Td(t)  =  f  g*s(t)  1/3 

T(t )  WsO)  9*s(tu) ) 


(7) 


The  success  of  big  bang  nucleosynthesis  (BBN)  and  the  structure  of  the  cosmic  microwave  backgromrd  (CMB)  power 
spectrum  place  tight  bounds  on  any  new  relativistic  degrees  of  freedom  in  the  dark  sector.  The  limit  on  the  effective 
number  of  light  neutrino  species  is  Nv  =  3.24  ±  1.2  at  the  95%  confidence  level  E3,  which  gives 

S*£4(^bbn)  =  \  x  2  x  (_ZV„  —  3)  <  2.52  (95%  confidence)  (8) 
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for  3  light  SM  neutrino  species  [35] .  The  5-year  WMAP  data  [23]  also  bounds  the  number  of  neutrino  species  by 
Nv  =  4.4  ±  1.5  at  the  65%  confidence  level,  and  the  7-year  WMAP  data  [30]  places  a  tighter  lower  limit  of  N„  >  2.7 
at  the  95%  confidence  level. 


A.  The  Boltzmann  Equation 


Let  us  assume  the  dark  matter  ip  is  a  stable  particle  that  annihilates  with  a  thermalized  annihilation  cross  section 
(crv).  The  general  Boltzmann  equation  governing  the  number  density  n  of  a  particle  of  mass  m  is 

n  +  3 Hn  +  (av)  (n2  —  n|q)  =  0  ,  (9) 

where  H  is  the  Hubble  parameter 


a 

and  riEQ  is  the  equilibrium  number  density 


47t3G<7* 


45 


(10) 


"EQ  “  pfy  jd3^E,Td  =  (i)  ■  <n> 

where  I\2  is  the  modified  Bessel  function  of  the  second  kind  of  order  two.  Generalizing  the  traditional  treatment,  we 
allow  for  the  possibility  that  the  mass  of  the  dark  matter  is  a  function  of  a  real  scalar  chameleon  field  </>  and 

denote  (^-dependent  masses  and  couplings  with  a  tilde. 


It  is  convenient  to  scale  out  the  effects  of  the  expansion  of  the  universe  by  defining 

Y  =  — 
s 


(12) 


{n.,j,{x)  and  Y(x)  are  taken  to  be  independent  of  <f>)  and  to  use  a  new  independent  variable,  related  to  the  cosmic  time 
t  through 


x(t)  = 


ttit 

W)  ’ 


(13) 


where  nriT  is  some  constant  mass  scale.  In  the  usual  derivation,  my  is  chosen  to  coincide  with  the  dark  matter  mass; 
however,  since  our  dark  matter  has  varying  mass,  we  use  this  constant  parameter  instead.  Defining 


b  = 


45 


47t3G  ijit 


(14) 
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allows  us  to  write 


dx  m,T 


rjtot 
7*  5 


dt  bx 

which  can  be  used  to  rewrite  the  Boltzmann  equation  for  the  dark  matter  as 

Y' Or)  +  -AY2  -  Ylo)  =  0  ■ 


Here  a  prime  denotes  a  derivative  with  respect  to  x ,  and 

2 7T2  g*s 


B  =  ( av ) 


45  .  / Otot 


&?7ly  , 


(15) 


(16) 


(17) 


which  may  depend  implicitly  on  cp  in  our  model  via  a  <f>  dependence  in  the  cross  section.  Note  that,  in  terms  of  these 
new  variables,  the  equilibrium  term  is 


*EQ  = 


45  g 


(2?r2)2g*s  \  mT 


X^)\CK,  ( 


£  mT 


(18) 


with  g  =  2  for  Dirac  dark  matter. 


It  remains,  at  this  level,  to  specify  Y(xi ),  the  initial  condition  for  Y.  We  consider  A  =  Y  —  Yeq,  the  departure  from 
equilibrium  |31j.  which  obeys 


B 


A'  =  -^eq  -  ^A(2Feq  +  A)  . 


(19) 


At  early  times  (1  <  x  <C  Xf),  Y  tracks  Yeq  extremely  closely  such  that  A  and  |A'|  are  small.  Note  that  in  the 
non- relativistic  approximation,  T  <C  fh^ (fa), 


yEQ  ^  x3/2 e-(x/()(m^(<t>)/mT)  , 

and  so  Y^q/Yeq  ~  —fh^{<j>) / (jnr  and  A'  ss  0.  Thus,  the  required  initial  condition  is 

x  v  t  \  ,  xfm^(fa) 

Y{xi)-Ym{Xi)+  , 

where  B((pfa  and  rh^(fa)  are  evaluated  at  the  initial  value  fa  =  faxfa. 


(20) 


(21) 


After  the  freeze-out  value  xf,  Y( x)  will  asymptotically  approach  a  constant  value  Yx.  The  energy  density  of  non- 
relativistic  dark  matter  today  is  then 


Po  —  fhy,  (xo)  (0o)Too 

2tt2 

= 'diifj{41o)Yoo  g*s,oT0  . 


(22) 


Having  generalized  the  usual  treatment  of  dark  matter  as  a  fluid  to  the  case  in  which  there  is  a  chameleon  field 
determining  the  dark  matter  properties,  we  now  turn  to  specific  examples  of  particle  physics  models  in  which  these 
phenomena  might  arise. 


III.  GAUGED  DARK  MATTER 


Consider  dark  matter  to  consist  of  a  Dirac  fermion  ip,  charged  under  a  dark  U(l)  gauge  group  with  gauge  boson  A^, 
and  a  dark  Higgs  field  $  that  spontaneously  breaks  the  U(l).  We  also  introduce  a  chameleonlike  field  <p  that  is  a 
real  scalar  field  with  properties  that  depend  on  the  dark  matter  energy  density.  The  chameleon  couples  to  the  other 
particles  in  the  dark  sector  by  entering  into  the  dark  matter  mass  m^(fa),  the  U(l)  coupling  /(</>),  and  other  couplings 
described  below.  We  consider  only  an  isolated  dark  sector  so  that  we  may  investigate  the  properties  of  this  simple 
model  without  the  complications  of  coupling  to  the  visible  sector. 
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A.  A  Toy  Model  for  Varying  Coupling 


As  a  first  step,  let  us  consider  the  QED  Lagrangian  with  a  real  scalar  held  <j>,  but  in  which  we  allow  the  coupling 
constant  e  to  vary  as  a  function  of  spacetime  (32j .  Specifically,  it  can  vary  as  a  function  of  <f> .  Let  us  write  the  new 
coupling  as  /(</>).  Thus, 


£qed<£  =  , 

z  4  J*(<p) 


(23) 


where  =  d^A^  —  d^A^.  Making  the  redefinition  A M  — >  /(0)Aa1,  we  obtain 

£qed^  =  -\d^cj)d^(t)  -V ((/>)  +  -  m^i\>  -  [^(/A")  -  3"  (/A")] 

z  4:jz  t  J 

Both  Lagrangians  are  equivalent,  but  now  the  gauge  transformation  reads 

/(0)Am  — »  f((j))A^  +  d^to 



e~lulifj 


+20; 


V’  • 


(24) 


(25) 

(26) 
(27) 


If  we  can  neglect  factors  of  (d^f/f)  compared  to  all  other  mass  scales  in  the  theory  (except  the  Planck  mass),  then 
the  Lagrangian  simplifies  to  the  approximately  gauge-invariant  form 


£qed0  ~  -V {<)>)-  -  m^ipip  -  f{<j>)^4^ 


(28) 


with  U(l)  current 


f{x)  =  . 


(29) 


B.  The  Cosmological  Equations  of  Motion 


We  now  include  gravity  and  a  complex  dark  Higgs  held  $  to  break  the  U(l)  symmetry  and  give  the  dark  gauge  held 
a  mass.  We  allow  for  a  varying  dark  matter  mass  by  using  the  effective  mass  parameter  my,  (<)>),  and  in  the  spirit  of 
effective  held  theory,  we  also  allow  all  couplings  [not  just  the  U(l)  coupling  /(</>)]  to  depend  on  cj>. 


Neglecting  factors  of  (d^f/f),  the  action  is  then 


S 


-  V(<t>)  -  (D^D^  -  H0(4>) 

-  -  my  (0)^0  -  Ay(0)($  +  4>t)^i/’  , 


(30) 


where  the  gauge  covariant  derivative  is  D ^  +  */(</>)AM.  The  equations  of  motion  for  the  helds  then  follow  as 

{iJp  -  -  Ay (</>){&  +  $*))  ip  =  0  (31) 

a<j)  -  v'((j>)  -  -  /,(^>)V’4'0  -  Ay(^)(4>  +  4>t)^V’  =  0  ,  (32) 


where  a  prime  denotes  differentiation  with  respect  to  <j>.  Let  us  assume  that  the  universe  is  dark-charge  symmetric,  so 
the  average  charge  current  density  is  negligible  compared  to  the  dark  matter  number  density  [see  (37)  below].  Thus, 
the  term  proportional  to  f'/f  should  be  small  compared  to  the  one  containing  to'/to,  given  that  f'/f  ~  to'/to  to 
within  a  few  orders  of  magnitude  -  a  condition  we  will  enforce  later.  We  may  write  this  last  equation  as 


Ucj)  -  V'{(j))  -  rh -  Ay (</>)($  +  $t)V>V’  «  0  • 


(33) 
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We  will  arrange  for  the  dark  Higgs  to  have  a  sufficiently  large  mass  that  its  perturbations  are  negligible  and  simply 
replace  $  by  ($)  in  the  equations  of  motion.  The  VEV  generates  an  additional  mass  term  for  ip,  but  we  can  redefine 
m^{(p)  to  absorb  this  term.  We  then  have 

(ilfi  —  rh!/,  ((/>))  ip  ~  0  (34a) 

0(p  —  V\cp)  —  m'^((p)ipip  «  0  .  (34b) 


We  calculate  the  energy-momentum  tensor  for  ip  by  varying  the  action  with  respect  to  the  metric.  Taking  care  to 
correctly  handle  the  nontrivial  metric  dependence  of  the  covariant  derivative  (33j.  we  have 


T$)  =  ['07(/xVI/)V>  -  (V^b^]  +  /(0)^7(m^)^  > 


^  2 

where  we  have  integrated  by  parts  and  used  the  field  equation  of  motion.  Taking  the  trace,  we  obtain 

i  r  -  .  -t—  i 
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—  __ 


flV 


2  L 

1 

2  L 


ip'f  ip  —  ip^}7  ip  +  f((p)ipj^ip 


y- 


ipi{ y  +  if((p)4)ip  -  ipi(f  -  if{(p)4)ip 

=  -rh^((p)ipip  , 


(35) 


(36) 


where,  again,  we  have  used  the  Dirac  equation  for  ip  and  ip  to  obtain  the  last  line.  If  we  model  the  dark  matter  as 
nonrelativistic  dust,  its  pressure  is  zero  and  so  the  trace  of  the  stress  tensor  is  approximately  given  by  —  p^.  Thus, 


Pp  =  fh^{(p)ipip  . 


As  a  final  step  in  this  section,  we  use  this  result  to  rewrite  the  cp  equation  of  motion  ( 34b )  as 

□0-Kff(0)  =  o, 

where  the  effective  potential  is 

Ves  =  V((p)  +  fh^((p)n^ 

=  V(<p)  +rh^{cp)Y(x)‘y^-g*s  (~^r)  ■ 


(37) 

(38) 


(39) 


IV.  CHAMELEON  BEHAVIOR 


With  a  complete  model  in  place,  we  now  turn  to  a  detailed  investigation  of  the  dynamics.  We  first  examine  the 
chameleon  field,  which  is  central  to  the  effect  we  seek.  Assuming  that  <p  is  homogeneous  and  isotropic,  so  that  we  can 
neglect  spatial  derivatives  in  □</>,  the  equation  of  motion  becomes 

cp  +  2>H(p  +  V'(<p)  +  rh'^((p)n^  =  0  .  (40) 

It  is  convenient  for  seeking  numerical  solutions  to  work  with  a  dimensionless  variable 


P  = 


rriT 


(41) 


and  to  use  x  as  our  independent  variable.  The  equation  of  motion  becomes 


rv//  \  2  ,  b2x2  dV 

P”{x)  +  -P'(  x)  + 


x  '  '  rn^glot  dip 


2  n2b2 


g*s 


dm 


(fi =Pmj' 


45x  glot  dcp 


Y{x)  =  0  . 


(42) 


(p—Prri'x' 


We  choose  the  initial  conditions  for  <p  to  begin  at  the  minimum  of  its  effective  potential  and  to  move  with  same 
initial  velocity  as  the  changing  minimum.  The  minimum  </>min  solves  the  equation  bffC^mm)  =  0,  so  one  of  the  initial 
conditions  for  this  equation  can  be  obtained  by  evaluating  this  expression  at  ay  ,  using  the  relevant  value  for  V(ay) 
from  (21).  Furthermore,  since  <pm -ln  is  a  function  of  x ,  the  initial  velocity  is  found  simply  by  taking  a  derivative  and 


using  the  Boltzmann  equation  to  obtain  the  relevant  value  for  Y'  (xP). 
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A.  Exponential  Potentials 


Our  goal  here  is  to  work  out  a  single  example  model  that  exhibits  the  effects  we  are  investigating,  while  at  the  same 
time  remaining  compatible  with  experimental  constraints.  For  simplicity  we  will  choose  exponential  functions,  which 
also  have  the  nice  feature  that  observables  approach  a  fixed  asymptotic  value  at  late  times. 

With  these  comments  in  mind,  we  therefore  choose  the  form  of  the  effective  potential  and  U(l)  coupling  to  be 

F(</>)  =  A4e~^/mi  (43a) 

m^((f>)  =m^  (l-  A2e~cl>/m2'sj  (43b) 

/(^)=e(l  +  A3e-^m3)“3,  (43c) 

where  A  and  to,/,  are  constants  with  dimensions  of  mass,  and  e  and  A2,  A3  >  0  are  dimensionless.  The  term  with 
A2  is  necessary  to  incorporate  the  properties  of  ip  into  the  equation  of  motion  for  cp.  The  possibility  for  A3  =  0 
(constant  gauge  coupling)  is  viable,  but  we  are  specifically  interested  in  increasing  the  cross  section  for  ip  as  the 
universe  expands.  We  choose  this  form  for  /  so  that  both  the  annihilation  and  scattering  cross  sections,  which  we 
calculate  below,  increase  with  time. 


The  largest  energies  of  the  particles  in  our  theory  are  of  order  m ^  for  nonrelativistic  dark  matter,  since  all  other 
particles  should  be  lighter  than  the  dark  matter  to  allow  for  annihilation.  We,  therefore,  require  m2,  m3  m ^  to 
suppress  higher-dimensional  operators  involving  derivatives  of  to ^  and  /  when  we  expand  the  action.  Additionally, 
we  need  to  1  >  A  to  suppress  higher-dimensional  operators  in  the  self-couplings  of  (p. 


The  effective  potential  in  ( 39 )  is  now 


VeS(<P)  =  A4e-^mi  +  m*(l  -  A2e-+/m2)Y(x)^g*s 


possessing  a  critical  point  at 


(44) 


TOito2  , 

- In 

m2  —  TOi 


toi  m^mj*  Y  \ 
m2  A4  a;3  J 


(45) 


which  is  real  and  finite.  In  order  to  generate  a  mass  for  the  excitations  of  cp,  we  require  this  critical  point  to  be  a 
minimum,  which  holds  for 


TO2  >  TOi  . 


The  minimum  moves  with  a  speed 


d^min 

dx 


TO1TO2 
TO  2  —  TOi 


ydF  3 

dx  x 


which  is  positive  (</>min  increases  with  x).  Finally,  we  identify  the  initial  conditions  for  </>: 


(p(Xi)  =  - 


TO1TO2 
m2  —  mi 
m\m2 


In  A 


TOi  TOy,TOjn  Y[xp) 


m2 


A4 


^(®»)  =  - 

dX  777-2  —  ml 



TO1TO2  / _3_  rh^  Teq(^)  +  / (4i?^TOr)\ 

to2  -  mi  \Xi  £ tot  YEq{xi)  +  (x^m^,) / {2B£,mT) ) 


In  order  to  ensure  rh^,  >  0,  we  require 


<p>  m2  ln(A2) 


(46) 

(47) 

(48) 

(49) 

(50) 


for  all  <p  relevant  for  our  calculation. 


B.  An  Attractor  Solution 


A  particularly  interesting  and  simple  possible  evolution  for  the  chameleon  field  is  for  it  to  begin  at  the  minimum 
of  the  effective  potential  and  then  to  adiabatically  track  this  minimum  as  it  evolves  cosmologically.  This  attractor 
solution  is  achieved  if  the  physical  mass  of  the  chameleon  satisfies 


^(f),  ph 


KWmin)  »  H  . 


(51) 


If  (51 1  holds  during  radiation  dominance,  when 


Hr  = 


(52) 


then  we  can  avoid  solving  the  coupled  differential  equations  ( 16 )  and  ( 40 )  and  simply  use  the  expression  for  <f>m in  for 
the  evolution  of  <p.  Similarly,  if  (51)  holds  during  matter  domination,  when 


HM  =  H0^y/2  ,  (53) 

then  we  can  easily  determine  the  value  of  <j>  today,  which  is  needed  to  calculate  the  values  of  the  ^-dependent 
parameters  today. 


Under  the  approximation  that  m2  mi, 


^0,ph 


27t2 

45  m\vri2 


1,2  Y'Rg'Js\-iR  . 


(54) 


It  follows  that  Hr  decreases  more  rapidly  than  TO«^ph  with  time,  whereas  during  matter  domination,  Hm  and  m^iPh 
have  the  same  x  dependence.  We  shall  verify  later  that  these  attractor  solutions  exist  by  numerically  solving  all  the 
relevant  equations  of  motion. 


V.  PARTICLE  PHYSICS  INTERACTIONS  AND  CONSTRAINTS 


In  the  adiabatic  regime  described  above,  we  now  have  all  the  ingredients  necessary  to  understand  the  cosmological 
evolutions  of  the  fields.  We  next  turn  to  the  particle  physics  phenomenology  of  the  model.  To  do  this,  we  rewrite  the 
action  ( 30 )  without  gravity  to  give  the  Lagrangian 


1 


c  »  -  -  V{<t>)  -  (DM$)T (£>"$)  -  VoW 

-  -  m^((^)#  -  /(^)V’4V’  -  (<£)($ 



(55) 


with  =  d,j,  +  if(4>)A 


M- 


A.  Breaking  the  Dark  U{  1)  Symmetry 


The  potential  of  the  dark  Higgs  field  $  is  chosen  so  that  this  field  acquires  a  vacuum  expectation  value  (VEV) 


(0|$(x)|0) 


V2  ‘ 


(56) 


Decomposing  <I>  into  two  real  scalar  fields  via 
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we  can  then  use  unitary  gauge  \ix)  =  0  to  rewrite  the  kinetic  term  for  $  as 

-  =  —d^hd^h  - \f2mv  +  h)2A»A^  . 


(58) 


Thus,  the  Goldstone  boson  is  eaten  to  give  the  dark  U{  1)  gauge  boson  A M  a  mass  Ma{4>)  =  The  Yukawa  term 

generates  a  contribution  to  the  mass  of  ip,  but  since  ip  already  has  a  Dirac  mass,  we  need  not  rely  on  the  dark  Higgs 
to  be  the  sole  source  of  the  ip  mass.  We,  therefore,  absorb  the  dark  Higgs  contribution  into  the  definition  of  to ^  and 
retain  the  freedom  to  choose  this  mass  scale  and  the  coupling  X^{(p)  separately. 


A  typical  choice  for  the  pure  dark  Higgs  potential  Vo(<p)  is 


W) 


\^h{4>) 


- y2 


which,  when  expanded  about  the  VEV,  yields 

V0(h)  =  \\h{(j>)\i2h2  +  yA h((p)\ih3  +  ^-X h{(p)h4  . 
4  4  16 

The  mass  of  the  physical  dark  Higgs  particle  h  is  therefore 


(59) 


(60) 


-  ,  Xh((p) 


and  we  see  that  the  masses  of  the  A  and  h  fields  are  then  related  by 

Ma{4>)  =  j2—  rhh(<p)  . 

V  Xh{(p) 


(61) 


(62) 


Since  the  relative  sizes  of  f(<p)  and  A h{<P)  are  unrestricted,  in  principle  the  relative  masses  of  A  and  h  are  not  fixed. 
However,  in  order  to  simplify  the  analysis,  we  will  impose  the  hierarchy  rhh{(f> )  >  2Ma{<P)  for  all  relevant  <p  so  that  h 
has  a  tree-level  decay  channel  to  A. 

Our  Lagrangian  at  this  stage  is  then 

£  =  -  ^d^cpd^cp  -  V{<p)  -  ^d^hd^h  -  *  A h(cp)y2h2  -  ^Xh(cp)vh3  -  ^Xh{(p)h4 
-  +  iip$ip  -  rh^(cp)ipip  -  f(cp)ip'yfJ‘Aflip 


i  r 
4  L 


2/2(0) 


h2A»A„  -  - 


1  r 


2  L 


2/(0)Ma(0)1  hA^A^  -  y/2)^,(if>)hipip  . 


(63) 


What  remains  is  to  incorporate  the  fact  that  </>  is  adiabatically  tracking  the  minimum  of  its  effective  potential.  To 
achieve  this,  we  expand  4>(x)  =  4>c(t)  +rj(x)  around  its  classical  value  and  recall  that  m2  and  m.3  are  sufficiently  large 
to  suppress  non-relevant  terms  of  0{rj)  or  higher.  The  Lagrangian  (63)  then  becomes 


£  =  -  \d^ri  -  ^hd»h  -  +  iPpU'  - 


V(^c)+2  V"(<Pc)r,2  +  0(V3) 


1  ‘ 
4  . 


2  L 


Xhi^  +  Oiv) 


1 


2  7  2 
M  h  — 


A^A^- 


-Xh{(pc)  +  0{r]) 


Mh3  - 


24 


f{<t>c)  +  o{ v)  tPi^A^ip 


-  [m^((pc)  +ifi^(4>c)ri  +  0(’i]2)]  ipip-  V2  X^{<j>c)  +  O(rj) 


4  l 


2f2(4>c)  +  0{rf) 


h2A^A„  -  l 


2f{<t>c)MA{<Pc)  +  0(rj) 


,  Xh  (<t>c)  +  0(ri) 

hipip 
hA»A„.  . 


(64) 
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FIG.  1:  Tree-level  ip  annihilation  diagrams.  The  massive  vector  boson  A  is  a  wavy  line,  and  the  scalar  h  is  a  dashed 
line.  Annihilations  to  A  +  A  and  h.  +  h  via  77-exchange  and  annihilations  to  final-state  77  particles  are  suppressed  by 

large-mass  factors. 


B.  The  Dark  Matter  Annihilation  Cross  Section 


Our  central  goal  is  to  understand  how  the  dependence  of  dark  matter  cross  sections  on  the  chameleon  field  changes 
the  standard  dark  matter  creation,  evolution,  and  detection  story.  To  this  end,  we  next  turn  to  the  calculation  of  the 
dark  matter  annihilation  cross  section.  The  relevant  Feynman  rules  can  be  found  in  the  Appendix. 


We  assume  that  the  dark  matter  is  the  heaviest  particle  in  the  dark  sector,  such  that  rh ^  fhh,  Ma-  Then,  the 
lowest  order,  tree-level  processes  for  2  — >  2  dark  matter  annihilation  are  shown  in  Fig.  [T]  and  their  amplitudes  are 


Mi  =  ie^e  2>v2 


f2(<pc)  -  k[) 7^  +  7AtAv,(pi  -  k'2)^v) 


V2f((pc)\,p((pc)MA((pc)Ah(pi  +P2)g , 


fiiy 


Ui 


(65) 


M2  =  iv  2 


2A|(0C)  (A^(pi  -  k[)  +  A^(pi  -  k'2)) 


+  —pz\^((pc)\h{(pc)v&-h{pi  +  P2) 


u  1 


(66) 


M3  =  ieVvv2 


V2f((pc)X^((pc)  (Av,(pi  -  k[)Y  +  7" A y,(pi  -  k[)) 

+  f2{<Pc)MA{(pc)%i^A  (Pi  +P2 ) 


U 1 


(67) 
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FIG.  2:  Ladder  diagrams  for  dark  matter  annihilation  (left)  and  scattering  (right). 


Working  in  the  center-of-mass  frame  and  in  the  nonrelativistic  limit  we  then  obtain 

(68) 

(69) 

(70) 

where  v  is  the  relative  velocity  of  the  incoming  particles. 

The  interaction  of  main  interest  is  between  the  dark  matter  and  the  gauge  boson  mediator.  The  dark  Higgs’s  primary 
role  is  to  break  the  U(l)  symmetry  to  give  the  mediator  a  mass,  and  most  of  its  particle  interactions  can  be  neglected. 
The  contribution  a2  is  p  wave  and  thus  subdominant  to  the  other  processes,  which  are  s  wave.  Also,  the  diagrams 
involving  exchanges  of  h  in  A4i  and  A42  do  not  significantly  contribute.  Thus,  the  dark  Higgs  has  the  opportunity 
to  influence  dark  matter  annihilations  only  via  03.  Let  us  insist  that  the  Yukawa  coupling  A ^  is  small  enough  (recall 
that  the  dark  matter  does  not  rely  on  this  coupling  to  obtain  a  mass)  such  that  <73  can  be  safely  ignored.  In  this  case 
only  cr  1  remains  and,  since  it  is  an  s-wave  cross  section,  it  is  a  simple  task  to  carry  out  the  thermal  averaging  required 
in  the  Boltzmann  equation.  Note,  however,  that  if  thermal  averaging  is  needed  (following  Ref.  [35] ) ,  we  must  use  the 
dark  sector  temperature  Td  in  the  expression 

(<™>  =  1  eqL  n2  oToATJ  f  ds  VsTdKi  ( (s  -  Am%)a{s)  .  (71) 

i.n^(Td))2^r  Jifhi  \TdJ 


<JlV 


a2v 


a3v 


/4(0c) 

167 

15A^(</>C) 

1287 Trh^((t>c)V 

m) A^e) 

87 rM|(</>c)  8ttv2 


C.  Corrections  to  the  Cross  Section 


We  are  interested  in  nonrelativistic  dark  matter,  for  which  the  relative  velocities  are  much  less  than  the  speed  of 
light.  It  is  well  known  that  for  sufficiently  low  velocities,  nonperturbative  effects  can  have  a  large  impact  on  the 
annihilation  and  scattering  cross  sections;  and  ladder  diagrams,  such  as  the  ones  shown  in  Fig.  [2j  must  be  included 
in  the  calculation. 


1.  The  Annihilation  Cross  Section 


In  the  case  of  annihilation,  performing  this  summation  is  equivalent  to  solving  the  Schrodinger  equation  in  quantum 
mechanical  scattering  theory  [36] .  This  yields  the  so-called  “Sommerfeld  enhancement”  (37]  of  the  annihilation  cross 
section  (for  detailed  reviews  in  the  context  of  dark  matter,  see,  for  example,  [551440] ) .  We  consider  the  annihilation 
cross  section  hq  for  a  pointlike  interaction  near  r  =  0  in  perturbative  field  theory.  For  small  velocities,  the  attractive 
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Yukawa  potential 


V(r) 


®e-MAr 

r 


(72) 


where  a  =  f2(4>c)/ 47T,  distorts  the  wave  function  at  the  origin  and  cannot  be  ignored.  Including  the  potential  will 
enhance  the  annihilation  cross  section  to  a  =  a^Sk  by  the  Sommerfeld  enhancement  factor  Sk ■  Let  us  define  the 
dimensionless  parameters 


tv 


e-A 


v 

a 

Ma 

CXTTl'ip 


(73) 

(74) 


where  v  is  the  velocity  of  each  annihilating  particle  in  the  center  of  mass  frame.  In  the  case  of  a  massless  gauge  boson 
with  a  Coulomb  potential,  it  is  possible  to  solve  the  Schrodinger  equation  analytically  to  obtain  the  Sommerfeld 
enhancement. 


For  a  massive  gauge  boson,  the  situation  is  more  complicated,  since  the  attractive  potential  has  a  finite  range  that 
limits  the  enhancement  from  being  arbitrarily  large  for  very  low  velocities.  In  the  regime  ca  <C  e,„,  we  recover  the 
Coulomb  case.  At  the  crossover  point  ev  ~  ca  (or  equivalently  m^v  ~  Ma),  the  de  Broglie  wavelength  of  the  dark 
matter  becomes  comparable  to  the  range  of  the  interaction.  At  lower  velocities  with  ca  e2,  the  Yukawa  potential 
cannot  be  ignored.  As  v  — »  0,  the  de  Broglie  wavelength  increases  to  a  value  larger  than  the  interaction  range,  and 
thus  the  enhancement  saturates  at 


en 


olMa 

rh ^ 


(75) 


Furthermore,  zero-energy  bound  states  may  form  for  certain  values  of  ca,  giving  resonance  regions  with  larger  en¬ 
hancements  ~  ca/c?,  until  they  are  cut  off  by  finite  width  effects.  In  the  early  universe,  freeze-out  typically  occurs 
at  velocities  Vf  ~  0.3,  so  that  ev  >  1  and  the  Sommerfeld  enhancement  can  be  ignored.  Note  that  there  are  no 
enhancements  for  ca  >  1. 


To  find  the  thermally  averaged  cross  section,  taking  into  account  the  Sommerfeld  enhancement,  we  integrate  Sk  using 
a  Maxwellian  distribution 


F(v)  =  ^±=v2e~^  , 

V6y/l T 

where  v  is  the  characteristic  velocity  of  the  astrophysical  system  of  interest.  Thus, 

(av)  =  (cru)s_waVe  (Sk) 


dv  F(v)Sk  ■ 


(76) 


(77) 

(78) 


For  the  purposes  of  this  paper,  we  choose  to  work  in  the  ca  >  1  regime.  This  has  two  consequences.  Practically,  the 
calculation  becomes  much  simpler,  since  we  need  not  worry  about  the  Sommerfeld  enhancement  at  all.  In  addition, 
by  deemphasizing  the  Sommerfeld  enhancement,  we  clarify  the  extent  to  which  the  novel  effects  developed  in  this 
paper  can  alone  increase  the  cross  section  over  time  in  areas  of  parameter  space  that  the  Sommerfeld  enhancement 
cannot  reach. 


2.  The  Scattering  Cross  Section 


To  find  the  scattering  cross  section,  we  can  use  nonrelativistic  quantum  mechanics  and  sum  over  partial  waves.  The 
total  cross  section  is 

A  _°°. 

0-  =  yy  J^(2 1  +  1)  sin2  Si  , 

1=0 


(79) 
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although  a  more  useful  quantity  to  compare  to  observational  constraints  is  the  transfer  cross  section 

/J 

dfl  (1  -  cos0)^ 

=  jyy  y, [(2/  +  1)  sin2  (5;  -  2(1  +  1)  sin  (5;  sin  Si+1  cos(<5;+i  -  Si)]  ,  (80) 

i 

which  controls  the  rate  at  which  energy  is  transferred  between  colliding  particles.  Following  (41],  analytic  estimates 
for  the  cross  section  are 


chr  — 


4-7T 

^Vrel 

4-7T 


(1  +  Lf 

(81) 

(1  +  L)  , 

(82) 

where  L  =  ^vTeibmax  is  the  largest  angular  momentum  needed  to  describe  the  interaction  between  two  particles  of 
reduced  mass  /i  =  m^,/ 2  that  travel  with  a  relative  velocity  ryei  and  maximum  relevant  impact  parameter  6max.  Note 
that  these  estimates  are  only  valid  for  L  >  1.  We  estimate  the  impact  parameter  by  solving 


1 

2 


A^rel 


f2/  4tt 


g-MAlmax 


(83) 


If  we  work  in  the  ca  >  1  regime  to  avoid  Sommerfeld  enhancements,  then  we  will  also  tend  to  avoid  enhancements  to 
the  scattering  cross  section  and  can  expect  to  be  working  in  the  Born  limit.  Simply  taking  the  nonrelativistic  limit  of 
the  perturbative  cross  section  gives 


/4(</>c)m^(^>c)  m^c) 

a  — - - - =  — - — i —  .  (84) 

8t TM\{<j>c)  8ttv4 

Assuming  that  dark  matter  self-interactions  are  not  needed  to  explain  the  structure  of  dwarf  galaxies  ED  ,  we  use  a 
conservative  bound  [32]  (see  also  IS  El) 

a/m y,  <  0.1  cm2/g  (85) 

for  characteristic  velocities  of  10  km/s.  As  we  mention  below,  it  would  not  be  difficult  to  find  parameters  that  violate 
this  bound. 

In  the  usual  treatment  of  dark  matter,  constraints  such  as  this  one,  obtained  from  present-day  observations,  can  be 
directly  applied  to  bounds  on  physics  at  freeze-out  or  before.  It  is  important  to  remember  here  that,  in  our  model, 
the  evolution  of  the  chameleon  field  means  that  such  a  connection  is  far  less  direct,  and  such  bounds  typically  do  not 
apply  in  the  early  universe. 


D.  Dark  Decays 


The  dark  Higgs  h  and  the  dark  gauge  boson  A  are  allowed  to  decay.  As  mentioned  earlier,  we  assume  m/j(0c)  > 
2Ma{4>c )  so  that  h  has  a  tree-level  decay  channel  to  A ,  as  shown  in  Fig.  [3)  Its  decay  width  is  then 

r  =  /2(^c)  mK&c)  /  4 

32tt  M\{(j>c)\l  rh\{<j>c) 

Although  the  A  particle  is  allowed  to  decay  to  rj  particles,  which  are  substantially  smaller  in  mass,  this  occurs  through 
a  1-fermion-loop  process,  as  shown  in  Fig.  |4j  The  amplitude  is  also  suppressed  by  two  factors  of  m2  from  the  0(r/2) 
term  in  the  expansion  of  {</>).  The  nonzero  amplitude  in  the  limit  of  frirj  <C  M a  ,  m^,  is 


1  - 


4Mi(^c)  ,  10MA(&) 


lMc) 


12- 



(86) 
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FIG.  3:  Tree-level  h  decay. 


FIG.  4:  1-loop  A  decay.  Only  the  first  diagram  is  nonzero. 


4i7r2/  , 

M  = - A - k'2  ■  e*(p)  Afh^MAC0]jp2 ,(p  -  k'2)2 

ma  L 

+(8to|  +  MA)B0[p2,rh^,rh^\  -  8m^B0[k'2,  m^,  rh^} 
where  Bq  and  C o  are  scalar  Passarino-Veltman  functions  pTTTl-l  l  (i| .  defined  via 


1 


1 


B()\p  jTTl  ,H~l  ]  •  2  /  ^  l  (12  i  2\[(1  |  \2  i  21 

mz  J  (lz  +  mz)[{L  +  p)z  +  mz\ 


Co[p2,{p-pi)2,p\,m1m.,m\  = [  d4l 

nrz  J 


1 


(l2  +  m2)[(l  +  p)2  +  m2][(l  +pi)2  +  m2] 
The  Co  integral  is  finite  and,  in  the  approximation  rh^  MA  mq ,  reduces  to 

1 


Co\p  ,{p  ~  k'2)  A  ,rh^,rh^,rh^}  «  ~^2 



(87) 

(88) 

(89) 

(90) 


The  Bq  integral  diverges,  so  we  cut  off  the  loop-momentum  integral  at  some  large  scale.  Using  m3  for  this  purpose, 
since  we  will  often  find  it  numerically  to  be  the  largest  mass-suppression  scale  in  our  theory,  we  have 


.B0[p2,?rty,77ty]  «  B0[k'2,  rh^,  m$\  «  2  In 

Putting  everything  together,  the  decay  width  of  A  is  then  given  by 
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f  m^,A-. 
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(91) 


(92) 


The  A  bosons  must  decay  efficiently  enough  not  to  contribute  significantly  to  the  energy  density  budget  today.  Though 
the  decaying  exponential  makes  meeting  this  criterion  more  difficult,  there  is  still  a  small  sample  of  parameter  space 
for  which  the  A  energy  density  does  not  pose  a  problem. 


VI.  NUMERICAL  SOLUTIONS 


While  we  have  described  a  number  of  ways  to  understand  the  evolution  of  the  fields  analytically,  including,  for  example, 
the  adiabatic  approximation  in  which  the  chameleon  tracks  the  minimum  of  its  effective  potential,  ultimately,  we  are 
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FIG.  5:  Scan  of  parameter  space.  Blue  points  indicate  sets  of  parameters  that  satisfy  all  constraints,  except  (for 
most  points)  for  having  a  negligible  A  energy  density.  Red  points  do  not  satisfy  the  scattering  cross  section  bound 
a/m t/,  <  0.1  cm2/g.  Green  points  do  not  satisfy  the  adiabatic  condition  in  (51 1  and  should  be  solved  with  the 

coupled  differential  equations. 


able  to  numerically  solve  the  relevant  equations  of  motion  completely.  To  do  so,  of  course,  we  must  make  sensible 
choices  for  our  parameters  to  satisfy  the  various  bounds  and  simplifying  inequalities  we  have  specified. 

We  need  to  implement  the  correct  relationship  between  the  dark  sector  temperature  Td  and  that  in  the  photon  sector 
T,  which  in  turn  requires  us  to  correctly  enumerate  the  massless  degrees  of  freedom  at  the  relevant  scales.  At  the 
unification  scale,  all  the  dark  particles  (ip,  A ,  h,  cp)  are  relativistic,  so  g^g(tu)  =  8.5.  Around  the  epoch  of  dark 
matter  freeze-out,  only  ip  is  nonrelativistic,  so  =  5-  Thus,  at  freeze-out,  £/  =  1.19  for  g*s(tf)  =  106.75  or 

=  0.56  for  g*s(tf)  —  10.75.  With  these  numbers,  the  bound  on  the  number  of  effective  neutrino  species  in  ([8])  is 
easily  satisfied. 

The  model  is  insensitive  to  Iff \{(p)  and  fhh(<p)  at  lowest  order.  We  choose  v  such  that  ipip  — >  AA  is  kinematically 
allowed  today,  while  ensuring  ca  >  1  and  a  <  1.  We  must  then  check  that  Ma(4> o)  satisfies  scattering  cross-section 
bounds.  It  is  simplest  to  assume  the  attractor  solution  for  (p  and  then  later  verify  that  it  is  in  fact  adhered  to.  The 
A  gauge  bosons  need  to  decay  away  before  BBN  so  that  their  energy  density  is  negligible.  Finally,  we  must  ensure 
that  the  evolution  ends  with  the  observed  density  of  dark  matter  today.  For  this  figure  we  use  the  bounds  from  the 
7-year  WMAP  data  ED],  assuming  a  ACDM  cosmology 

fWh2  =  —  =  0.1109  ±  0.0056  .  (93) 

PcO 


Given  these  constraints,  we  numerically  solve  the  Boltzmann  equation  and  show  a  sample  of  parameter  space  in  Fig.  [5j 
resulting  from  a  random,  uniform  scan  over  £  [0.1,500]  GeV;  mi  £  [105,107]  GeV;  m2  €  [5  x  105,  5  x  10s]  GeV; 
m3  €  [5  x  10s,  5  x  108]  GeV;  A  €  [10, 103]  GeV;  A2  €  [0.1, 9.9];  A3  £  [0.1,10];  and  e  £  [0.01,  V^tt]-  The  upper-left 
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FIG.  6:  The  number  of  A  decays  per  particle  between  freeze-out  and  BBN.  Points  above  the  horizontal  line  at  1 
indicate  that  all  A  particles  should  have  decayed  and  thus  do  not  contribute  significantly  to  the  energy  budget  of  the 

Universe. 


panel  shows  the  coupling  parameter  e  vs  the  dark  matter  mass  parameter  to,/,.  The  upper-right  panel  shows  the 
annihilation  cross  section  crv  vs  the  scattering  cross  section  a/rh ,/,,  both  evaluated  at  Xq  today.  The  bottom  panels 
show  the  boost  in  annihilation  cross  section  from  freeze-out  to  today  and  the  scattering  cross  section  today  vs  the 
mass  parameter  to,/,.  Again,  there  is  flexibility  when  choosing  v  without  affecting  the  evolution  of  <fi  and  Y  at  lowest 
order,  so  it  is  possible  to  obtain  valid  models  for  a  scaled  value  of  a/m ,/,.  Here,  we  show  the  largest  possible  scattering 
cross  sections,  while  staying  within  the  bound  6  a  >  1.  As  demonstrated  in  Fig.  [Gj  only  a  small  portion  of  the  sampled 
parameter  space  fulfills  the  requirement  that  the  A  gauge  boson  energy  density  is  negligible  by  the  time  of  BBN. 
While  finding  a  set  of  parameters  that  satisfies  all  constraints  is  certainly  possible,  the  effect  of  having  very  large 
increases  in  the  annihilation  cross  section  does  not  seem  to  be  a  general  feature  of  the  model. 

As  a  concrete  example,  we  show  a  specific  model  with  the  parameter  choices:  to ,/,  =  123  GeV;  tot  =  to,/,;  mi  =  38  TeV; 
m.2  =  500  TeV;  m3  =  500  TeV;  A  =  18  GeV;  A2  =  0.6;  A3  =  9.2;  e  =  0.96;  and  v  =  10  GeV.  This  comprises  an 
optimistic  set  of  parameter  choices,  that  satisfies  all  our  bounds  and  provides  a  large  change  of  order  ~  106  in  the 
annihilation  cross  section  over  the  history  of  the  Universe.  Our  choice  for  the  value  of  v  gives  €a  =  1-07  today,  and 
we  can  ignore  the  Sommerfeld  enhancement.  Larger  values  of  v  work  equally  well;  increasing  v  increases  £a  ~  v  and 
decreases  a/fh^  ~  v~4.  The  dark  matter  relic  density  is  =  0.1097,  within  a  standard  deviation  of  the  observed 
value.  The  scattering  cross  section  today  is  4.9  x  10-4  cm2/g,  well  below  the  conservative  limit  in  ([85]) .  We  must 
also  check  that  these  parameters  satisfy  the  assumptions  we  have  made  in  writing  down  the  model.  For  example,  we 

neglected  terms  with  d^f / /,  and  here  we  note  that  ///  ~  10-9  —  10-6  GeV,  which  is  much  smaller  than  other  mass 
terms  in  the  perturbative  expansion.  The  adiabatic  approximation  is  satisfied  with  H/m^p h  ~  10-11  throughout  the 
evolution  of  (/>.  Finally,  we  use  the  decay  width  of  the  A  particles  to  determine  that  they  have  decayed  away  in  the 
time  from  freeze-out  to  BBN,  so  they  do  not  contribute  to  the  energy  budget  we  observe  from  the  CMB. 

The  results  for  the  evolution  of  <j>,  V,  the  dark  matter  mass  to,/,,  and  the  coupling  /  as  a  function  of  T  =  mj* / x 
are  shown  in  Fig.  [7]  We  also  show  the  annihilation  and  scattering  cross  sections  in  Fig.  [8j  The  scattering  cross 
section  quickly  approaches  its  asymptotic  value  by  the  time  of  dark  matter  freeze-out,  while  the  annihilation  cross 
section  still  grows  orders  of  magnitude  from  freeze-out  to  now.  This  difference  is  due  to  the  scattering  cross  section, 
a/fh^  oc  771^,/v4,  and  the  annihilation  cross  section,  av  oc  f4/m^,  depending  differently  on  <p  via  and  /.  We 

choose  the  form  of  /  to  force  the  annihilation  cross  section  to  grow  more  slowly,  whereas  the  scattering  cross  section 
has  no  such  term  countering  its  growth.  With  these  particular  choice  of  parameters,  the  scattering  cross  section  is 
too  small  to  have  interesting  astrophysically  observable  consequences. 

As  shown  in  Fig.  [5]  there  are  other  choices  of  parameters  that  will  still  give  a  boost  to  the  annihilation  cross  section 
while  yielding  a  larger  scattering  cross  section  to  match  observational  bounds  [ill];  however,  again,  most  of  the 
plotted  parameter  space  is  restricted  from  the  A  energy  density  requirement.  One  option  for  increasing  the  viable 
parameter  space  is  to  relax  the  requirement  that  €a  >  1  and  to  work  in  the  regime  of  Sommerfeld  enhancements; 
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FIG.  7:  The  evolution  of  0  (upper  left),  Y  (upper  right),  fh^,  (lower  left),  and  /  (lower  right)  as  a  function  of  T  in 
GeV  for  the  model  with  m ^  =  123  GeV.  The  red,  dotted  line  indicates  the  approximate  dark  matter  freeze-out 

temperature. 


our  model  would  still  provide  significant  increases  in  the  cross  sections,  and  Sommerfeld  enhancements  would  serve 
to  further  increase  the  boosts.  Another  clear  option  is  to  open  an  alternative  decay  channel  for  A. 


VII.  CONCLUSIONS 


In  this  paper  we  have  investigated  the  possibility  that  the  properties  of  dark  matter  depend  crucially  on  the  dynamics 
of  a  chameleon  field  -  a  scalar  field  whose  cosmological  evolution  depends  not  only  on  its  bare  potential  but  also  on  the 
local  density  of  other  matter  (such  as  dark  matter  itself)  in  the  Universe.  We  have  shown  that  such  a  coupling  allows 
the  annihilation  cross  section  (for  example)  of  the  dark  matter  particles  to  change  by  several  orders  of  magnitude 
between  freeze-out  and  today,  while  remaining  consistent  with  all  observational  constraints.  We  have  presented 
a  general  formalism  to  describe  how  this  might  happen  and  have  provided  a  specific  particle  physics  example,  in 
which  all  relevant  quantities  can  be  calculated.  While  there  are  significant  observational  and  theoretical  constraints 
on  models  of  this  type,  it  is  nevertheless  possible  for  the  cross  section  to  evolve  in  such  a  way  that  there  may  be 
interesting  implications  for  the  detection  of  dark  matter  and  for  its  dynamical  effects  on  late-universe  astrophysics. 

There  are,  of  course,  other  possible  complications  to  this  idea  that  are  beyond  the  scope  of  the  current  paper  but  that 
provide  interesting  avenues  for  future  study.  One  natural  step  is  to  couple  our  model  directly  to  the  Standard  Model. 
One  way  to  achieve  this  is  to  directly  add  the  dark  U(l)  to  the  current  SM  gauge  group  [32].  Another  possibility  is 
to  couple  to  the  Standard  Model  through  1/(1)  kinetic  mixing  [48!,i49j.  This  extension  of  our  model  should  be  able 
to  easily  accommodate  the  relevant  particle  physics  constraints  [5DH55].  while  easily  allowing  for  decays  of  the  dark 
gauge  boson  to  Standard  Model  particles  well  before  BBN.  The  dark  matter  annihilations  would  still  be  dominated 
by  the  channel  00  — ►  A  A,  since  annihilation  to  Standard  Model  particles  would  be  suppressed  by  the  small  coupling 
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FIG.  8:  The  evolution  of  the  annihilation  cross  section  (av)  (left)  and  the  scattering  cross  section  a/rh ^  (right)  as  a 
function  of  T  in  GeV  with  =  123  GeV.  The  red,  dotted  line  indicates  the  approximate  dark  matter  freeze-out 

temperature. 


parameter  for  the  U(l)  mixing.  However,  it  is  a  more  delicate  issue  to  decide  what  a  natural  route  would  be  to  couple 
the  visible  and  dark  scalar  sectors,  particularly  with  regards  to  coupling  the  chameleon  to  normal  matter. 

Finally,  we  did  not  attempt  a  careful  analysis  of  the  effect  of  late-universe  inhomogeneities  on  the  chameleon  field  or 
the  dark  matter  properties  on  which  it  depends.  In  the  specific  models  we  considered,  it  seems  as  if  such  effects  would 
be  small,  but  a  more  careful  examination  is  warranted. 
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Appendix:  Feynman  Rules 


The  Feynman  rules  are  shown  in 


Fig-®  All  of  these  diagrams  have  higher-order  corrections  that  involve  77  particles. 


FIG.  9:  Feynman  rules  for  h  (dashed  line),  ip  (solid  line),  A M  (wavy  line),  and  77  (dotted  line).  We  include  the 
Yukawa  interaction  with  ip  and  ? 7,  which  is  relevant  for  the  1-loop  A-decay  amplitude  in  (87),  but  other  77-interaction 
vertices  are  not  shown.  All  parameters  labeled  by  a  tilde  are  evaluated  at  (j)c. 
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